Testing excision techniques for dynamical 3D black hole evolutions 
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I. INTRODUCTION 

A successful numerical model of a binary black hole col- 
lision will require the refinement of a number of special- 
ized techniques. A particular problem of black hole sim- 
ulations is handling the strong fields in the neighborhood 
of the physical singularity. The method of "singularity 
excision" promises to solve, or rather avoid, this diffi- 
\ culty. Black hole excision was first discussed by Thorn- 
burg based on a suggestion by Unruh Q. Within 
the event horizon of a spacetime all light-cones corre- 
sponding to the flow of physical information in the fields 
point "inwards" . Thus a boundary within the event hori- 
zon forms a one way membrane, with information flow- 
ing only into the hole. Physical data on such a boundary 
should not influence the external spacetime. In a numer- 
ical simulation, this implies that errors on the bound- 
ary, for instance due to extrapolations from the exte- 
rior, should not be expected to propagate to the exte- 
rior spacetime. However, because of the complicated and 
nonlinear nature of the equations being solved, complex 
relationships link the finite differenced evolution equa- 
tions to artificial boundary conditions. It is often the 
case that the full implications of using a particular tech- 
nique are poorly understood, and can only be judged by 
some amount of experimentation. 

The excision idea has a long history of implementa- 
tions in numerical relativity codes in ID, 2D and 3D. 
The basic concept was developed into a set of techniques 
for practical application in ID in 0, where causal dif- 
ferencing techniques and dynamical lapse and shift con- 
ditions designed for excision were introduced and shown 
to increase stability and accuracy of dynamic BH space- 
times (coupled to self-gravitating scalar fields) by orders 



of magnitude. This work was extended both to 3D BH 
evolutions |^ and to a large class of shift conditions in f^l 
such as the geometrically motivated distance freezing, 
area freezing, expansion freezing and minimal distortion 
shift. These ideas were later a pplie d to more complex 
spacetimes, e.g. H H IHIll El 111 El , typically using 
the Arnowitt-Deser-Misner (ADM) evolution system and 
analytic shift conditions. 

With the switch to the more stable BSSN evolution 
system [H El El El 

Il9j | and the accompanying im- 
provements in stability of 3D codes, the excision idea was 
revisited by Alcubierre and Briigmann in the form of an 
algorithm which they called "Simple Excision" 0, ■ 
The idea was to attempt an implementation which re- 
moved some complications that had hindered previous ef- 
forts, namely causal differencing within the horizon, and 
complicated interpolations onto an irregular boundary. 
Instead, the excision region was fixed as a cubical region 
with faces along constant Cartesian coordinate planes, 
simplifying the determination of a normal, and a sim- 
ple source-term-copy boundary condition replaced true 
causal differencing. The result proved to be remarkably 
successful at evolving distorted single black hole space- 
times allowing accurate evolutions of unlimited length in 
time. 

More recently advances in excision in 3D Cartesian co- 
ordinates have been reported, using different tec hniq ues, 
for scalar fields in fixed black hole backgrounds |22, 123 
(see also |24| ) and for evolutions of single black hole 
spacetimes |l3l l25l l2^ . In general, the excision problem 
becomes much better controlled when using coordinates 
that are adapted to the excision surface, and furthermore 
when using hyperbolic formulations and making use of 
the characteristic information about the propagation of 
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all degrees of freedom. Spherical excision is used in spec- 
tral codes j23,|2^429,J but has also been proposed for finite 
difference codes using multi- patch techniques [snLlsiIls^ . 

A basic premise of the excision method is that the ex- 
cision boundary data should not propagate to the exte- 
rior spacetime. In practice, it is important to verify that 
this is in fact the case for a given numerical implemen- 
tation, and to quantify any errors that might be intro- 
duced. In this paper, we propose to compare long term 
black hole evolutions using singularity excision to corre- 
sponding evolutions carried out over a full grid, but with 
otherwise identical initial data and evolution methods. A 
useful excision method should yield essentially identical 
evolutions, and thus identical physical measurements, in 
the region of the spacetime outside of the horizon. 

The main candidates for evolution without excision are 
methods using singularity avoiding slicings, where the 
gauge condition on the lapse slows evolution of the hy- 
persurface in regions where a physical black hole singu- 
larity is approached. In fact, singularity avoiding slic- 
ing conditions (for example, maximal slicing or algebraic 
lapse conditions of the "1-1- log" family) were the method 
of choice in most simulations prior to the development 
of excision techniques. Without appropriate shift condi- 
tions, slice stretching arises due to the differential infall 
of coordinates, severely limiting the evolution time. It 
has only very recently been shown that an appropriately 
chosen shift vector can alleviate these problems, so that 
they occur at a much later time (or even not at all for 
a fixed resolution) and the final black hole settles to its 
expected stationary state 13 ■ 

For this reason, comparisons between evolutions with 
and without excision are rarely done. Even worse, in 
most cases where excision has been successfully used, 
very crude measures, such as "time to crash" , or some 
overall measure of the Hamiltonian constraint have been 
used as a measure of success. With the exception of a 
small number of results [1TII2H . it has generally not been 
verified that, for example, correct gravitational wave- 
forms can be extracted with excision, a point which is 
crucial for upcoming gravitational wave observations. As 
these waves will typically be very small perturbations (of 
order 10~^ or less) of the metric [s^ [Sa, very small 
errors generated by the excision techniques could easily 
swamp the signals being extracted. Furthermore, even 
if the excision properly preserves the causal structure, 
allowing no physical signals to propagate out from in- 
side the horizon, gauge effects may well propagate to the 
outside, influencing the solution. This could show up 
not only in metric functions but also in so-called gauge- 
invariant waveforms, since such waveforms are only in- 
variant under infinitesimal gauge transformations given 
some assumptions about the background coordinate sys- 
tem. A large gauge wave propagating through the space- 
time may very well be seen in waves extracted in this 
way. 

In this paper, we propose to use both waveforms, ex- 
tracted far from the horizon, and detailed information 



extracted from apparent horizons just outside the exci- 
sion region, as important test quantities for spacetimes 
evolved both with and without excision. Such tests would 
have been difficult or impossible in the past because, un- 
til recently, 3D BH simulations without excision became 
rapidly inaccurate and crashed after evolution times of 
only t « 30 — AOMadm, where Madm is the ADM mass. 
However, the recent development of powerful shift con- 
ditions has to some extent cured the problem of slice 
stretching associated with singularity avoiding slicing, 
making it possible to evolve certain classes of distorted 
and colliding BH spacetimes much longer (for thousands 
of Madm) and much more accurately than ever before, 
without the need for excision ^33] . This makes it possible 
for the first time to carry out systematic, long term stud- 
ies of the effects of the excision technique on the evolved 
BH spacetimes. 

Following Refs. [M IM El (sl we use distorted black 
hole initial data sets, whose waveforms can be indepen- 
dently and reliably computed, as testbeds for the exci- 
sion techniques we use. Those papers showed that even 
modes with energies of order 10^'^ Madm could be very 
accurately extracted in a full 3D numerical evolution, as 
confirmed by comparisons with purely perturbative evo- 
lutions. We will use some of the same testbeds proposed 
there in this paper. 

We also apply the ideas to the study of BH collisions, 
using Brill-Lindquist-type initial data which, as showed 
in Ref. (33j, can be evolved indefinitely without resort- 
ing to excision. Although Misner data have been much 
more extensively studied [3^ 1^ and would be an 
excellent test-bed, we use Brill-Lindquist here because of 
their close relation to the Brill-Lindquist family of "punc- 
ture" data often applied in generating data for binary 
BHs in quasi-circular orbit 0, EE IM EE E^ • These 
data sets possess apparent horizons in the initial slice, 
allowing excision regions to be defined, but have also 
been evolved for unlimited times without excision (af- 
ter merger). As such, these evolutions provide an ideal 
test-bed for excision methods. If the excision method is 
valid and the exterior spacetime is protected from errors 
on the excision boundary by the event horizon, then sim- 
ulations with excision should accurately reproduce the 
non-excised results. Because the testbeds we use here 
are highly dynamic, distorted and colliding black holes, 
evolved in general coordinate systems in 3D Cartesian co- 
ordinates, without resort to special conditions or tricks, 
successful tests of excision techniques against indepen- 
dently computed physics results would provide a certain 
amount of confidence that the techniques can be applied 
to more complex systems, such as orbiting black holes, 
where the detailed physics is more critical to get right. 

In the next section we briefly outline the methods used 
in our black hole simulations. There follows a description 
of the excision region and boundary conditions which 
we are applying. Finally, we present numerical results 
in which distorted and binary black hole evolutions us- 
ing excision are compared directly with those performed 
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without excision. B. Lego Excision 



II. MODEL AND METHODS 

A. Initial Data and Evolution Methods 

The details of the formulation which we have imple- 
mented numerically are outlined in 33J . Briefly, we make 
use of the "BSSN" evolution system [H [H E^. The 
evolution variables are a conformal factor, the confor- 
mally decomposed metric, the trace-free extrinsic curva- 
ture tensor and its trace. Additionally, the contracted 
conformal Christoffel symbols := 7-''^r*^ are evolved 
as independent variables. These variables can be re- 
garded as gauge source functions, and in our evolutions 
are used to fix the shift vector through e.g. the condi- 
tion QiF* = 0, which results in an elliptic shift condition 
analogous to minimal distortion 47]. In practice, instead 
of applying this elliptic shift condition directly, we use a 
hyperbolic "Gamma driver" version as described in Eq. 
(46) of m, 

df(3^=Fdtr-{rj-^)dtf3\ (2.1) 

The parameter 77 is a driver term which controls the 
growth of the shift and is discussed in more detail later, 
and -F is a function which can be used to condition its 
overall shape. 

The black hole initial data we use comes in two classes. 
The first class is the "Brill wave plus black hole" family 
of distorted BHs, constructed topologically as wormholes 
isometrically connecting two asymptotically flat space- 
times. These data sets have proved a rich and powerful 
system for developing many techniques in numerical rela- 
tivity and for studying the physics of many aspects of dis- 
torted BHs 21, .35, 36, 37, 48, 49, 50, 51, 52, 53, 54, 55]. 

The second class is Brill-Lindquist binary BH initial 
data in isotropic coordinates '56'] . The topology of these 
data sets is that of an asymptotically flat spacetime con- 
nected through wormholes to two disconnected asymp- 
totically flat ends. These data sets are time symmetric 
and set the stage for BH head-on collisions. 

As discussed in Is^ , for spacetimes in which there 
is enough symmetry (such as the axisymmetric head-on 
collision) , it is possible to carry out expansions of the evo- 
lution of the metric quantities at the "punctures" (the 
points in TZ^ at the "center" of each BH representing 
the compactified asymptotically fiat regions on the other 
side of the wormholes), thus ensuring that the evolution 
variables remain regular there. By careful consideration 
of the evolution at the punctures, it is therefore possi- 
ble to carry out long-term evolutions with the punctures 
included on the grid, i.e. without excision. 



The excision boundary condition which we apply is 
a variation of the "Simple Excision" methods described 
m The important difference is that instead of excis- 
ing a cubical region, we excise an irregular region. The 
size and shape of the excision region is determined dy- 
namically by the location of the apparent horizon. On 
a Cartesian grid, the resulting surface is a blocky quasi- 
spherical region, prompting the name "lego excision" for 
this method. 

The update of all interior points on the grid (including 
those within the apparent horizon) is carried out using 
centered differences. That is, no causal differencing is 
applied. The exceptions are advection terms on the shift 
(terms of the form f3^di ), for which a second order up- 
wind in the shift direction is applied. 

For points on the excision boundary itself, a centered 
update scheme cannot be used for lack of data to one side. 
Instead, data at these points are updated using a time 
derivative term which is copied from one point out. This 
has the effect of allowing a certain amount of evolution 
to take place on the boundary, while also allowing the 
boundary to "settle down" as the system settles to what 
should be a stationary end state, exemplified by dtcj) = 
for the fields 67] . The neighboring point from which 
the time derivative is copied is determined by a normal 
to the surface at the point. The time derivative from the 
point nearest to the normal is copied to the boundary 
point. 

For the models in question, we assume that the exci- 
sion region neither shrinks nor moves across the grid. In 
other words, it is not possible for grid points to emerge 
from the excision region: once excised, a point remains 
excised. While this may at first seem like a strong limi- 
tation, in fact even with these restrictions the technique 
is adequate for many binary BH inspiral studies. Under 
reasonable gauge choices it should be expected that the 
BH horizons individually should either remain of an ap- 
proximately fixed coordinate size, or grow slightly, as is 
the case for the final post-merger BH. Further, for binary 
BH inspirals, co-rotating coordinates reduce the amount 
of dynamics on the grid to the dynamical time-scale of the 
problem in question, and either remove or greatly reduce 
the orbital motion of the individual bodies. Co-rotating 
coordinate systems have been proposed (see, e.g. 58]), 
and implemented for cases of BH binaries |59ll60l l61] (for 
neutron stars see e.g. H^)- For the puncture data which 
we have used here, the punctures are fixed to the grid. As 
such, problems involving movement of the excision region 
do not arise in these cases, almost by construction. 
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FIG. 1: The absolute value of the Hamiltonian constraint 
\H\ for a distorted black hole simulation without excision, 
evaluated along the y-axis at T = 9.38Madm {x = z — 0). 
The high resolution Hamiltonian is scaled such that the two 
curves should coincide for second order convergence. 



FIG. 2: Absolute value of the Hamiltonian constraint \H\ 
along the j/-axis for the same distorted black hole simulation 
shown in Fig. but this time using excision. The Hamilto- 
nian constraint was computed at T = 9.38Madm ■ 



III. APPLICATIONS 
A. Results for Distorted Black Hole 

As a first application we study a simulation of a dis- 
torted BH which has been studied extensively in previ- 
ous work m m m El S S 113. The BH is ini- 
tially distorted by an even parity Brill wave and, for a 
low enough amplitude wave, during the evolution the BH 
rings at quasi-normal frequencies before settling down to 
a Schwarzschild BH. The parameters used here are, in 
the notation of [Ulillillei, Qo = 0.1, % = 0, cr = 1. 
For this simulation the computational domain extends to 
±34.80 (in coordinate units). For these initial data pa- 
rameters, Madm = 1.92, which puts the outer boundary 
at ±18.13Madm- We also use two different resolutions, 
the smallest one with a 288 x 288 x 144 sized uniform grid 
and a coordinate resolution of 0.24, and the largest with 
a 384 X 384 x 192 sized uniform grid and a resolution of 
0.18. An explicit reflection symmetry about the z = Q 
plane is used. 

We evolve this system both with and without exci- 
sion, using identical evolution parameters. For the gauge 
we use a Gamma driver shift condition as specified in 
Eq. ()2.1|l . and 1-t-log slicing. In the simulation using ex- 
cision, the lego excision region was located at 80% of the 
apparent horizon radius, with a minimum buffer size of 
at least 5 grid points between the excision region and 
the apparent horizon surface. We find the horizon u sing; 
methods described in 0,113, and implemented as in [ggl^ 
At higher resolution the size of the buffer was adjusted 
such that the excision region remained at approximately 
the same position. 

In Fig. n and Fig. El we show convergence plots for 
the absolute value of the Hamiltonian constraint along 
the y-axis at time T = 9.38Madm, for the runs without 



and with excision, respectively. In both plots the higher 
resolution case is scaled by a factor of (0.24/0.18)^ so 
that, in case of second order convergence, the two curves 
should coincide. Note that, in both cases, we have sec- 
ond order convergence in the regions that are causally 
disconnected from the outer boundaries. In Fig.|21points 
inside the excision region have been removed from the 
plot. Note that outside the excision region the curve is 
very similar to that of Fig. The noise visible in both 
plots at y > TMadm is due to numerical details of the 
initial data constraint solver. 

In Fig. 13 we plot the difference between the Hamilto- 
nian constraint between runs with and without excision 
at T = 9.38M^£)Af at the two different resolutions, with 
the high resolution case again scaled as in the previous 
plots. Only the points outside the excision region are 
shown in the plot. The figure shows that there are in- 
deed small differences between the two cases, but that 
these differences converge away at second order outside 
of the excision region. 

As an indicator for analyzing excision effects we use 
the ratio Cr of the circumference along a polar direction 
in the j;2-plane, to the circumference in the equatorial xy- 
plane, of the apparent horizon, as a function of time. As 
these measurements are made very close to the excision 
region, and as for the low amplitude Brill wave studied 
here Cr remains within about 1% of the Schwarzschild 
value of unity, we regard this as a very sensitive indicator 
for these tests. In Fig. ^ these ratios are plotted for 
evolutions with excision (lines) and without (dots) at two 
resolutions. The differences are very small in comparison 
to the size of the physical wave, and much smaller with 
increased of resolution. Note that in this plot the non- 
excision curves do not extend over the same time as the 
excision data. The earlier termination of the non-excision 
runs in Fig. 01 was due to a hardware problem and was 
not caused by problems in the code. 

In Fig. [SI we show the evolution of the Zerilli wave 
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FIG. 3; Differences in the Hamiltonian constraints between 
the runs with and without excision — Hi\o — cxc| at T — 

9.38Madm for the distorted black hole. The difference be- 
tween the two cases converges away at second order. 



FIG. 5: Zerilli function {£ = 2,m = mode) calculated 
at r = 13.5Madm ■ The lines show the excision data and the 
points show the results of the simulations without excision. 
The inset shows the absolute value of the difference in the 
waveforms which decreases as the resolution is increased. 
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FIG. 4: Ratio Cr of the polar and equatorial circumference 
of the apparent horizon for two different resolutions of the 
distorted black hole simulation. The lines denote the excision 
data and the dots show the results of runs without excision. 
The inset shows the absolute value of the difference between 
the two cases. 



function ■0 (£ = 2, m = mode). The figure shows 
excision data as lines and the results from simulations 
without excision as points. The absolute value of the 
difference between the waveforms is also shown in the 
inset. This difference is much smaller than the Zerilli 
function itself and decreases with resolution. 

We also studied the effects of the size of the excision 
box for this model. One can expect that problems will 
occur if the excision box is made too large and in par- 
ticular if it reaches outside of the black hole. In Fig. El 
we show the influence of differently sized excision boxes 
on the waveform. For these simulations the black hole 
is nearly spherical (see Fig. ^ with a mean coordinate 
radius of initially about 1 which then grows to about 2 
during the first ^M^^jjm of evolution. For comparison 
we have also used a cubical excision box which remained 



fixed throughout the evolution. The resolution used in 
these simulations was 0.18. We first plot the waveform 
of the lego excision simulation used above. Next we 
show waveforms for cubical excision regions where half 
the length of one side of the cube is 0.54, 0.90, 1.26 and 
1.98, respectively. For the case with the smallest excision 
box the waveform remains essentially unaffected. How- 
ever, once the excision region becomes larger and even 
extends outside of the initial apparent horizon a large 
effect on the waveform becomes visible. This demon- 
strates that excision must be done inside the horizon as 
one would expect, ft also indicates that with techniques 
implemented here, a buffer zone of two to three points 
in the initial hole can be sufficient to accurately extract 
gravitational waves. 



B. Results for Head-on Collision 

As a second example we consider a binary BH head- 
on collision for Brill-Lindquist type data, starting from a 
close coordinate separation between punctures of 2.303. 
For this data set each individual BH has a mass parame- 
ter of Mi = 0.5 for a total ADM mass of Madm = 1 5^ 
The initial data uses a "fish-eye" coordinate transforma- 
tion to push the boundaries further out ^i,|6Q|- We have 
performed runs using octant symmetry with central res- 
olutions of 0.128, 0.064, and 0.032, on cubical grids with 
96'^, 192"^ and 384'^ grid points respectively. Using the 
"fish-eye" parameters a = 3, s = 1.2 and tq = 5.5, in 
the notation of [sj, places the coordinate boundaries at 
12.288 and the physical boundaries at 25.862. 

The evolution without excision is identical to the one 
shown in Figs. 11 and 12 of Ref. [s^. However, run- 
ning the excision case with exactly the same parameters 
turned out to be impossible. The reason is that with 
the choice of damping coefficient rj used in that reference 
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FIG. 6: Zerilli Function i/) for tiie distorted blacf: hole calcu- 
lated at r = 15.3Madm , using a resolution dx = 0.18 and 
different cubical excision box sizes. For the first two evo- 
lutions, the excised region is within the horizon and show 
essentially identical results for both the cubical and "lego" 
excision region. The other evolutions have at least a portion 
of the excision box extending outside the horizon, and thus do 
not represent true outflow boundaries. They show dramatic 
differences from the previous waveforms. 

for the Gamma driver shift condition, the lapse and shift 
turned out to be too dynamic at the excision boundary, 
causing the simulation to crash very early. We have found 
that this effect can be controlled by using a smooth spa- 
tially varying 7/-parameter. This was implemented using 
the conformal factor from the initial data 

with Ml and M2 the mass parameters of the black holes 
(in this case Mi = M2 — 0.5) and ri and r2 the coordi- 
nate distances to the two punctures. A smoothly varying 
77-parameter can then be constructed (this construction 
is of course not unique) as 

_ ??punc ~ Voo ,„ „N 

'^"''P""^" C(1 + (V|/-1)2)' ^■^■^> 

where jypunc is the value of rj near the punctures, rjoc 
is the value of r] at infinity, and C is a parameter that 
controls the width of the transition region. Here we used 
??punc = 5.6, rjac = 2.8 and C = 1. 

The resulting waveforms for the excision runs, using 
the above construction for r] and the three different res- 
olutions, can be seen in Fig. These waveforms are 
second order convergent as can be seen from Fig. |S1 
where the difference between the low and medium res- 
olution waveforms and four times the difference between 
the medium and high resolution waveforms are plotted 
together. 

At all three resolutions we can see some differences in 
the waveforms between the runs with and without ex- 
cisions. There are basically two reasons for this. The 
first is of course the presence of the excision boundary 



FIG. 7: Zerilli function il) {£ — 2, m = mode) calculated 
at r = IA.SMadm for the head-on collision using excision at 
resolutions 0.128 (low), 0.064 (med) and 0.032 (high). 
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FIG. 8: Scaled difference of the Zerilli function -0 (^ = 2, 
m — mode) calculated at r = 14.8Madm for the head-on 
collision using excision showing second order convergence. 



and the second is the slightly different gauge choices that 
were needed in each case in order to obtain long enough 
evolutions. However, as can be seen from Fig. |21 these 
differences converge away to second order with resolu- 
tion (except for a small gauge effect which is visible in 
the initial part of the waveform). 

For the above runs, we used an excision region which 
was 80% of the apparent horizon radius with minimum 
buffer size of (5,10,20) grid points for the three different 
resolutions. This was in order to ensure that we had 
a buffer zone large enough suppress inaccuracies at the 
excision boundary from propagating out of the hole, as 
well as to ensure that the apparent horizon can always 
be located, as this requires points on both sides of the 
surface. 

We note that the runs performed with excision are 
more sensitive to gauge parameters than the evolutions 
without excision. This is essentially due to the fact that 
the shift can be used to control the rate of growth of 
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FIG. 9: Scaled difference in waveforms from evolutions with 
and without excision for the three different resolutions used 
in the head-on collision. 

the horizon, and in some cases can even draw it inwards. 
For the excised runs, this means that the horizon can 
be brought too close to the excision boundary, allowing 
errors to "leak out". The shift parameters used in this 
paper allow a slig htly faster growth of the horizon than 
those used in 1331. 



IV. CONCLUSIONS 

We have proposed a sensitive test-bed for 3D black 
hole excision techniques in numerical relativity, whereby 
simulations with and without excision are carried out on 
identical initial data sets. A recently developed family 
of powerful gauge conditions now allows long evolutions 
of some black hole spacetimes without excision, making 
such a test-bed practical for the first time. Applying 
these ideas to our own methods, we demonstrated that 
our excision procedure is not only effective in improving 
the length of time that black holes can be evolved, but 
also that it preserves sensitive details of the physics that 
can be extracted. We carried out both distorted and bi- 
nary black hole head-on collision simulations, both with 
and without excision. We used both the ratio of apparent 
horizon circumferences and the very sensitive indicator of 
extracted waveforms to compare physically important de- 
tails of the simulations. In particular, the waveforms rep- 
resent very small signals buried in the metric functions. 
The results are essentially identical, indicating that the 



excision boundary condition has preserved the accuracy 
of the runs. 

The evolutions were carried out for cases of axisyni- 
metric initial data, where for singularity avoiding slicing 
with an appropriate shift condition the evolution vari- 
ables could be suitably controlled at the locations of the 
punctures. For more general situations it can be more 
difficult to maintain a regular evolution. In such situa- 
tions, even for the puncture data discussed above, it is 
expected that excision would be necessary for a long-term 
evolution. However, the tests shown here have given us 
confidence that the excision methods which we are using 
are robust for reasonably complicated situations, and do 
not adversely affect the accuracy of evolutions. 

Runs without excision form an important benchmark 
against which the effects of excision techniques can be 
analyzed. Evolutions of black hole spacetimes involve a 
patchwork of experimental techniques, all of which are 
needed to obtain long evolution times, but each of which 
carries with it complications that can affect the accuracy 
and stability of the simulation, both in isolation and in 
interaction with the entire system. Given this fact, the 
models considered above can form a particularly impor- 
tant test of the effects of applying an excision boundary 
condition. 

The initial data sets we used here are readily available 
to anyone through the Cactus framework. This will make 
it possible for other groups to apply the same test-bed to 
their own evolution codes. 
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